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Signatures of Quantum Chaos in Wave Functions Structure for 

Multi-well 2D Potentials 
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We propose a new approach to investigation of quantum manifestations of classical stochasticity (QMCS) 
in wave functions structure, which can be realized in potentials with two and more local minima. The main 
advantage of the proposed approach is the possibility to detect QMCS in comparison not different wave func- 
tions, but different parts of the same wave function. Efficiency of the approach is demonstrated for two 
potentials: surface quadrupole oscillations (QO) and lower umbillic catastrophe (UC) D5. 
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Energy spectra and eigenfunctions of classically non- 
integrable systems represent the main object of search 
for QMCS PmEI. It should be pointed out that in 
■analysis of QMCS in the energy spectra the principal 
role was given to statistical characteristics, i.e. quan- 
tum chaos was treated as property of a group of states. 
In contrast, the choice of a stationary wave function as 
'a basic object of investigation, relates quantum chaos 
.to an individual state. Usual procedure of search for 
QMCS in wave function implies investigation of distinc- 
tion in its structure below and above the classical energy 
i0f transition to chaos (or other parameters of regularity- 
'chaos transition) . Such procedure meets difficulties con- 
.nected with necessity to separate QMCS from modifica- 
'tions of wave functions structure due to trivial changes 
in its quantum numbers. Up to present time correlations 
'between peculiarities of the classical motion and struc- 
ture of wave functions were studied mostly for billiard- 
type systems 01 E] . For Hamiltonian systems with 
non-zero potential energy QMCS were studied either for 
model wave functions [Jj or for potential energy surfaces 
(PES) with simple geometry [H]. Till now there is prac- 
tically no information on wave functions structure for 
generic Hamiltonian systems, including multi-well po- 
tentials. Such systems allow existence of the mixed state 
(MS): different (regular or chaotic) classical regimes co- 
exist in different local minima at fixed energy ll(J| . 
Aim of our work is to show, that such systems repre- 
sent optimal object for investigation of QMCS in wave 
functions structure. Wave functions of MS allow to find 
QMCS in comparison not different eigenfunctions, but 
different parts of the same wave function, situated in 
different regions of configuration space (corresponding 

^'e-mail: cherkaskiy@kipt.kharkov.ua 



Figure 1: The level lines of the QO potential IBbll for 
= 13 (a), = 18 (b), and for the UC D5 Q with 
a = 2 (c). 



to different local minima of the potential). 
Let us demonstrate this possibility for MS, generated 
by the deformation potential of surface QO of atomic 
nuclei and lower UC D5 IT . It can be shown |13 |. 
that using only the transformation properties of the in- 
teraction, the QO potential takes the form 



C/Qo(ao,a2) 



m.n 



)'"aQ(6a2 



«o)" (1) 



where Oq and 02 are internal coordinates of nuclear sur- 
face undergoing the QO 

R{0,ip) = Ro{l+aoY2.oi0,ip)+a2[Y2.2{0,^)+Y2,-2{0,ip)]} 

(2) 

Since in the construction of ^ only transformation 
properties of interaction play role, this expression de- 
scribes potential energy of surface QO of a charged 
liquid drop of any nature (for example, a metal clus- 
ter ,14 ), containing specific character of the interaction 
only in the coefficients Cmn- Restricting ourselves with 
the terms of fourth order in the deformation and assum- 
ing equality of masses for the two independent direc- 
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tions, we get the following Can-symmetric Hamiltonian 
jj^P^^±PiL + UQo{x,y;a,b,c) (3) 



2m 



where 



UQoix,y;a,b,c) = 
f (x2 +y')+b _ Lx') +c{x' + y'f (4) 
X ^ao,y = V2a2,a = 2Cio,6 = 3Coi,c = C20 

Let us introduce the dimensionless variables 

{x,y) = lo{x,y), {px,Py) ^po{px,Py),E ^ eoE (5a) 



(5b) 
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In the variables {x, y) (further we will drop the bar line) 
the Hamiltonian J^l has the form 



H='^-^p^ + UQo{x,y;W) 
2m 



(6a) 



UQo{x,y;W) 



2W 



(x'^ +y'^) + xy^ 



'■ + {x' + y') 



y 

(6b) 

Hamiltonian Ht)a|) and corresponding equations of mo- 
tion depend only on = b'^jiac), which is the unique 
dimensionless parameter, that can be constructed from 
a, b and c, and it completely determines the PES 
[Fig^a),(b)]. Region < < 16 includes potentials 
with only one critical point — minimum in the origin 
[Fig^a)] , corresponding to spherically symmetric equi- 
librium shape of the nucleus (or liquid charged drop). 
For W > IQ the PES has seven critical points: four min- 
ima (one central and three peripheral) and three sad- 
dles, separating the peripheral minima from the central 
one [FigH^b)]. In this Letter we consider in details the 
case = 18, when the potential Ij6b|l has four minima 
with the same value Emin = and the saddle energies 
Es = 1/20736. It was shown [Sj, that the critical energy 
of transition to chaos E^r has different values for differ- 
ent minima: Ecr = Egjl for the central minimum and 
E'er — Eg for the peripheral ones. It means, that for 
Es/2 < E < Es regular and chaotic trajectories coex- 
ist and are separated not in phase, but in configuration 
space, resulting in the phenomenon of MS. 
MS is a common case for multi-well potentials. Ac- 
cording to catastrophe theory, a wide class of multi-well 
2D polynomial potentials can be generated by germs of 
lower UC of types D5 , of , D7 from the Thom catas- 
trophes list, affected by certain perturbation To 



demonstrate the proposed approach we consider, apart 
from the QO potential Ij6bp . a lower UC I?5, described 



by germ x'^ /4 + y^x with perturbation bx^ — ay^ . This 
potential has only two local minima and three saddles 
[Fig^c)], and therefore it is the simplest potential, 
where MS is observed. Under the Maxwell condition 
b = the energies of all the saddles are the same, 
and energies of all the local minima too. We will con- 
sider the case a = 2. = 1 



UDrAx,y) = '— 



xy 



2y^ - x^ 



(7) 



when all Emin = — 1 and all Eg ~ 0, and MS is observed 
in the energy region — 1/2 < E]\fs < 0. 
Calculation of quasiclassical part of the spectrum for 
systems with multi-well PES requires appropriate nu- 
merical methods. Matrix diagonalization (MD) method 
is attractive for the Hamiltonians with eigenfunctions, 
that do not differ too much from the basis functions. 
However this numerical procedure becomes less attrac- 
tive (or even not efficient at all) at the transition to PES 
of complicated topology (multi-well potentials). In par- 
ticular, the diagonalization of the QO Hamiltonian (|6afl 
with > 16 in the harmonic oscillator basis requires 
so large number of the basis functions, that goes beyond 
the limits of the computation power. In this case the 
attractive alternative to the MD may become the spec- 
tral method (SM) in the form, proposed by Feit et al. 

m 

Numerical solution of the stationary Schrodinger equa- 
tion 

h"^ 1 

-Y'^ + U{x,y) ^pn{x,y) ^ Enipn(x,y) (8) 

by the SM requires computation of the correlation func- 
tion 

P{t)^ J dxdy^-^ix,y)4'ix,y,t) (9) 

where ip{x,y,t) represents a numerical solution of 
the corresponding time-dependent Schrodinger equation 
with an arbitrary initial condition ipo{x, y) = ip{x, y,t = 
0). The solution tp{x,y,t) can be accurately generated 
with the help of the split operator method 



^(x,y,i + At) = 
e^^^e-'^*^(-^^)e'^^V(x,2/,i) + 0(At3) 



(10) 



where exp(i/i^AiV^/4)^(x, y, t) is evaluated with the 
help of the band-limited Fourier series representation 

N/2 N/2 

H^,y,t)= E V'™„(i)e^('"^+"^' 

m=-N /2+I n=-N /2+1 

(11) 
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Figure 2: Determination of two close energy levels _E„ 
and E„+i in the QO potential (16 bt with W — 18 and 
different number of time steps M = 2*,^: = 9, 10. ..18: 
(a) - absolute value of P^''\E) in log scale, (b) - val- 
ues of E^!'^ and -E^^i, (c) - consecutive corrections 



e: 
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where TV is the number of grid points along a grid line 
and Lq is the grid length. 

From the other hand, the solution 'ijj{x,y,t) can be 
expressed as a linear superposition of eigenfunctions 
il)n{x,y) of ijHIl 

ifj{x,y,t) = '^anipn{x,y) exp{~iEnt/h) (12) 

n 

We assume no degenerate states in the decomposition 
H12|) . which can always be achieved by certain choice 
of the initial condition tljo{x,y). Using (|12(l in 0, wc 
obtain 

P{t) = J2 eM-^Ent/h) (13) 

n 

The Fourier transform of (|13|) 



1 

P{E) = i y dte'^ P{t)wit) = ^ Kl^SriE - E„) 



where 



4 MX. 
e h -1 
■ ET 



ET 

e ^ h 
-ET 



4/ ET 

1 _|_ e'*~fi~ 



(Tr-2-) 



(14) 



(15) 



and w(t) — 1 — cos(27rt/T) is the Hanning window func- 
tion. The plot for P{E) displays a set of sharp local 
maxima a.t E = En [FiglS^a)] , where iJ„ are the energy 
eigenvalues of Once the eigenvalues are known, the 
corresponding eigenfunctions can be computed by nu- 
merically evaluating the integrals 



i'n{x,y) = ^ / dtij{x,y,t)w{t)e'''^' (16) 



This procedure is very efficient, when implemented with 
the help of the fast Fourier transform algorithm (FFT), 



and very accurate, since the spatial derivatives are ap- 
proximated to iVth order in Ax — Ay = Lq/N. Lq 
and N must be chosen large enough, and ^Q(x,y) suf- 
ficiently fast decaying with its Fourier components, in 
order to assure that tp{x, y, t) is negligible on the grid 
boundaries both in coordinate and reciprocal spaces. 
The sampling interval At limits the spectral bandwidth 
AEmax = Trh/At of a function that can be repre- 
sented by a Fourier series determined by sampled val- 
ues. Therefore At should be chosen small enough to 
accommodate the necessary number of energy levels, or 
^E„iax > AUmax^ wherc AUjnax is the maximum ex- 
cursion of the potential. Since the potential (|6b|) is un- 
bounded, it is necessary to put an appropriate cutoff in 
order to apply the SM. If the wave function tp{x,y,t) 
is generated over a total time T = MAt, the mini- 
mum separation in energy levels that can be resolved 
is AEmin = nh/T, which also provides an estimate of 
the accuracy with which individual eigenvalues can be 
determined from the numerically computed P{E) (|14(l 
without the aid of lineshape fitting techniques. Such 
techniques, however, may improve the eigenvalues accu- 
racy by roughly two orders of magnitude. Fig[21 shows 
the typical shape of P{E) [Fig[2Ia)] and the correspond- 
ing two close energy levels [Fig|2Ib)], calculated for in- 
creasing M = 2'' , k = 9,10, ...18, using the single-line fit. 
For A; = 9, 10 the levels look as one, at fc > 11 they are 
resolved, and for fc > 13 their values change much less 
than the level spacing. Figl^c) shows the differences 
between the eigenvalues, calculated with given M = 2*^ 
and the more accurate, calculated with AI = 2*^+^ — 
the procedure converges fast indeed, displaying the sat- 
uration in accuracy. 

In applying the SM, a preponderant fraction of the com- 
puter time goes to the generation of the 'ip{x,y,t), and 
most of that time in turn is invested in FFT computa- 
tions, so the overall calculation time for our 2D prob- 
lems scales as MN'^lnN. We performed our calcula- 
tions with TV — 512 and M — 65536, which allowed 
to obtain for reasonable time about 10^ eigenfunctions 
with sufficient accuracy and high coordinate resolution 
to allow their detailed analysis. The scaled Planck's 
constant h is an arbitrary parameter and is chosen to 
obtain the desired number of energy levels. 
Let us now discuss the results obtained for the poten- 
tials Uqo and J/d^. In the former case for Vl^ < 16 
the only possibility to detect the QMCS in the wave 
function structure is to look how it changes as energy 
grows. As was shown in for 4 < W < 16 energy 
region of chaotic motion is bounded from both sides: 
Ecri < Ec < Ecr2, which means the regularity-chaos- 
regularity (R-C-R) transition at energies Ecri and Ecr2 
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R, C Rj 

Figure 3: The R-C-R transition in the QO potential Ij6b|l 
with W = 13. 



respectively. We will distinguish three energy regions: 
low-energy regular Ri{E < Ecri), chaotic C{Ecri < 
E < Ecr2) and high-energy regular R2{E > £'cr2)- 
For the considered case = 13 the R-C-R transi- 
tion is observed at critical energies Ecri — 8 • 10^^ and 
Ecr2 — 8.4 • 10~^ respectively. FigO shows the changes 
in the structure of level lines of IV'nP [FigOJa)] and the 
corresponding nodal domains picture [Fig|3|[b)], clearly 
correlating with the character of the classical motion, 
displayed in the corresponding Poincare surfaces of sec- 
tion [FigEIc)]. 

The considered possibility corresponds to the traditional 
approach in search of QMCS in the wave function struc- 
ture. Existence of the MS, at > 16 for the QO po- 
tential [FiglUa)] or in the UC potential [FigEfb)], 
opens a new possibility. Comparing the structure of 
the eigenfunction in central and peripheral minima of 
the QO potential [Fig^fc)], or in left and right minima 
of the UC Dc, potential [Fig^Jd)], it is evident that the 
nodal structures of the regular part and the chaotic part 
of the eigenfunction are clearly different: 

i) within the classically allowed region the nodal do- 
mains of the regular part of the wave function form a 
well recognizable checkerboard-like pattern T ; nothing 
similar can be observed for the chaotic part; 

ii) the nodal lines of the regular part exhibit crossings or 
very tiny quasicrossings; in the chaotic part the nodal 




Figure 4: The MS in the QO potential (a),(c) and in 
the Ds UC (b),(d): (a),(b) - Poincare surfaces of section, 
(c),(d) - nodal domains of the eigenfunctions. 

lines quasicrossings have significantly larger avoidance 
ranges; 

iii) while crossing the classical turning line U{x^ y) = 
En, the nodal lines structure of the regular part im- 
mediately switches to the straight nodal lines, going to 
infinity, which makes the turning point line itself easily 
locatable in the nodal domains structure; in the chaotic 
part an intermediate region exists around the turning 
line, where some of the nodal lines pinch-off, making 
transition to the classically forbidden region more grad- 
uate and not so manifesting in the nodal structure. 
In conclusion we remark that the Hamiltonian system 
with multi-well PES represents a realistic model, de- 
scribing the dynamics of transitions between different 
equilibrium states, including such important cases, as 
chemical reactions and nuclear fission. Existence of the 
MS must essentially determine dynamics of the phys- 
ical processes in such systems, for instance the quan- 
tum assistance tunnelling. We demonstrated a possi- 
bility to observe the QMCS in an individual quantum 
mechanical state — an eigenfunction of the MS. Fur- 
ther analysis implies investigation of the eigenfunction 
amplitude distribution, nodal lines quasicrossings avoid- 
ance range distribution, and the wave packets dynamics. 
Another interesting perspective is to study relevance of 
the Berry-Robnik formula for the energy levels spacings 
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distribution in the MS and to investigate the nodal do- 
mains and the nodal lines statistics. 
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